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Abstract 

111 an effort to disc'over tlie causes for disagreeiiient 
between previous 2-D computations and nominally 
2-D experiment for flow over the element McDon- 
nell Douglas 30P-^t0N airfoil configuration at high 
lift, a combined experimental/C'FD investigation is 
des(T‘ibed. The experiment explores several difleriTit 
side-wall boundary layer control venting patterns, 
documents venting mass flow rates, and looks at <’or- 
ner surface flow patterns. The ex|)erimental angle 
of attack at maximum lift is found to be sensitive to 
the side wall venting pattern: a particular pattern 
increases t he angle of attack at. maximum lift by at 
least 2®. A signihcant amount of spanwise ])ressure 
variation is present at angles of attack near maxi- 
mum lift. A CFD study using TD structured-grid 
com])utations, wliich iindudes the modeling of side- 
wall venting, is employed to investigate' .‘TD effects 
on the flow. Side- wall suedion strength is found to 
affect the angle at which maximum lift is predicted. 
Maximum lift in the (’FD is shown to be litnited 
by the growt h of an off-body corner flow vortex and 
('onsec|uent increase in spanwise pressure' variation 
and dc'crease in edrculation. The TD e'om]nitations 
with and without wall venting predict similar trends 
to experiment at low angles of att ae'k, but either stall 
too early or else overprexlict lift levels near maximum 
lift by as much as 5%. Unstructured-grid c'omputa- 
tions demonstrate that mounting brackets lower the 
lift levels near maximum lift conditions. 
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1 Introduction 

A large number of (TT) studies have been (*on- 
d noted for multi-elemenl airfoil configurat ions over 
the last decade. Some of these pertained to exper- 
iments conducted in the NASA Langley low turbu- 
ieiue pressure tunnel (LTFd ) for McDonnell Dou- 
glas Telement configurations such as the ;f()F-!U)N. 
See. for example, refs. 15. Tlie LLF I tests were 
nominal 1>^ 2-D, obtained with the use of a side- wall 
i)oimdary layer control system that a]>j)lied suction 
near the model via venting througl) jiorous plates.'' 
d he (TD studies employed 2-D <'alculat ions wlu'u 
('om paring with LTFd data. 

Although CFD could claim some success in the 
prediction of the experimental multi-element flow 
field, certain key predictions were in error.' Most 
notably, (Tl) consistently predieted both the max- 
imum lift and the angle of attack at whic'h it otx'urs 
to be higlu'r tlian experiment (e.g., at n — 25° as 
opposed to — 21°). Although turbulenc'e mod- 
els were originally sus])ected as a potential ('ause for 
this diserepancy, additional ('FD studies with more 
advaiu'cd models*-'^*^'^ indi('ated that the differeiu'es 
between turbulenee models tended to be relati\ely 
small for these (’ases. 

A recent ('FD control surface effectiveness study 
by Jiang*' demonstrated that even nominally 2-D 
wind tunnel experiments with side- wall vejiting ('an 
have significant 5-D effeed.s that are too large to ig- 
nore in ('FD modeling. In an effort to determine 
if 3-D effects are the root cause of the disagreement 
l)etween ( 'FD and ex]>erimenls for the 30F-30N near 
maximum lift , an investigation using a two- pronged 
approach was undertaken. 1 he first effort was an 
experimetital investigation that focused on the ef- 
fects of side- wall venting and venting pattern on re- 
sults near maximum lift, and included attempts to 
deduce flow features near the side walls using sur- 
face oil. Also, the ex]:>eriment measured venting 
paramet('rs. siuJi as mass flow rate, that were not 
measured in earlier ex[)eriments on this configura- 
tion. I'he second effort involved modeling the 3-D 
configuration with ( 'FD, including side- wall venting. 


1 

American Institute of Aeronautics and Astronautics 



Sensitivity to venting suction levels was investigated 
and the flow field features near the side wall were 
examined. Additionally, CFD using unstructured 
grids was used to explore the effe(’ts of the mounting 
brackets on the .*1-1) flow field. This ]>a[)er summa- 
rizes the results from this two- pronged study. 

2 Description of the Experiment 

The is a McDonnell Douglas Aerospace 

(MDA) ^belemeiit configuration (designated LB54G 
in MDA’s nomenclature). The model consists of the 
LSI 2 slat, the WTOBB wing box assembly (main), 
and the F22 flap. The slat is at an angle of at- 
tack of —.40®, with a gap of 2.959? c and overliaiig 
of —2.5% c. The flaj) is at an angle of attack of 30^, 
with a gaj) of 1 .27% r and overhang of 0.25% c. The 
configuration has a stowed chord length c of 22 in,, 
and the LTPl' tunnel width is 3G in. For all results 
to b(' showii below, th<‘ wing coordinates have been 
nondimensionalized by c. In this coordinate system, 
in the de]:>loyed ]>osition the leading edges of each 
element are located in the following positions: slat 
at r/r — — Q.Q85L main at r/c = 0.0438, flap at 
,r/c — 0.8715. riie 3()P-30N ('onfiguration has been 
tested in the past at several different Reynolds num- 
bers ranging from 5 million to 10 million, although 
the current ex [>eri mental study only used Re 9 
million, d'he nominal Maxd] number in the current 
investigation is AJ — 0.2. 

In an attempt to achieve the most 2-D flow field 
possible, ]>orous [)lates were located at the side-walls 
of the LTPd near the wing.^^^ The porous plates 
were <'onnected to a venting chamber that vented 
out to the atmosphere through a remotely-operated 
ball valve. 4 he default experimental procedure was 
as follows: during any .set of giv(uj runs. th<^ ball 
valvf' was adjusted (with the wing at an angle of at- 
ta(‘k of 10°) to minimize the spanwise pressure vari- 
ation, as determined from several rows of spanwise 
pressure taps on the model. Then, at other angles 
of attack during the same series of runs, the valve 
remained open th<' same amount. In the (’urrent ex- 
periment , t he effect of varying the venting mass flow 
rate at each angle of attack near maximum lift was 
al.so investigated. A no- venting configuration (all 
porous plates ('overed with tape) as well as 4 dif- 
ferent porous plate <*onfigurations were tested; these 
are summarized in Table 1. 

Sk('tches of the porous regions used \u the experi- 
ment are depicted in Fig. 1. The porous region over 
the main element was not varied. It was approx- 
imately U.5 inch high and began about 5 in. back 
from the element s leading edge on its upper sur- 
face and ext(UKled downstream until it ('onnected in 


1’able I. Summary of venting configurations in the 
currenl experiment 


Con fig. 

main 

venting 

flaj) 

venting 

strip in front 
of slat 

1 

JiO 

no 

no 

2 

ve^s 

yes, widejiing 

no 

3 

yes 

yes, widening 

yes 

4 

yes 

yes, constant 

no 

5 

yes 

yes, ('oust ant 

.ves 


a continuous fashion with the j)orous region oveT 
the flap. Two configurations over the flap were 
tested. The fiist (“widening") had the region over 
the flap widening from 0.5 inch to approximately 
1. 5 inch near the trailing edge. The second (“con- 
stant'’) employed the default pattern used in j)ast 
tests (unpublished, but documented in LTPT engi- 
neer test notes), which kept the width over tln^ flap 
constant at 0.5 inch. When employed, the porous 
region in front of the slat was a large semi-chrcular 
strip of width I inch, whose downstream inner radius 
tou('hed t he slat leading edge. The porous plate had 
0.0025 inch diameter holes spaced 0.191 inch apart, 
for a porosity level of approximately 3 *4%.. In the 
experiment at Re = 9 million, the LTPT operating 
total pressure was approximately 54 psi. For most 
of the runs, unless otherwise noted, the Ap bet ween 
the vent chamber and the tunnel (freestream) was 
approximately 8 psi. 

3 Numerical Method 

The compressible CFD codes used in the cur- 
rent investigation were CFL3D,‘‘ a structured-grid 
upwind finite-volume method, and FI N3D,‘“ 
an unstructured-grid upwind finite- volume method. 
The turbulence model employed by both codes was 
the Spalart-Allmaras (SA) model, version In 

CFL3D, transition was set on each of the elements 
at the following fixed locations for all runs: x/c = 
—0.026 1 and —0.0847 on the slat lower and upper 
surface, x/c = 0.0999 and 0.0082 on the main lower 
and upper surface, and x/c — i.l243 and 0.92 1 4 on 
the flap lower and upper surface. These locations 
are representative of locations measured by Bert(4- 
riid,^'" but were not varied with angle of attack for 
this study. In ( TL3D. transition was achiev^ed by ze- 
roing out the turbulence production term in the re- 
gion of each grid zone where laminar flow is desired. 
All FUN3D runs assumed fully-turbulent flow. 

In simulating wall venting (used for the CFL3D 
runs only), a linear form of th<' Darcy pressure- 
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velority law was iisod:"' 

r«: = — - — (Pc - /^r). (I) 

where pc is the vent chamber pressure and a is 
an "effective" geometric porosity factor. Tims, the 
wall venting suction velocity was assumed to be di- 
rectly proportional to the difference l)etween the 
vent chamber pressure (taken in this case to be con- 
stant) and the local pressure at the wall. However, 
there is still some uncertainty in the a])plication of 
this law to the current case of the wing in the LTPT. 
Idle factor a in Darcy's law is somewhat ad-hoc. It 
is ]>roportional to the actual ])orosity level but is 
considerably larger in magnitude^'' and is usually 
determined f>y ('omparison with measured wall ve- 
lo(dt ies. In this study a range of values for rr from U - 
O.d in ste])s of 0. 1 were used to det ermine their effect . 
Resulting mass flow t hrough t he vent was compared 
with experimentally measured levels to ascertain the 
approximate correspondence bet ween the (TD and 
the experiment. 

In the CFD runs, the nondimensional vent cham- 
ber pressure was assumed to be p/ (p^ ) — 0.()2, 

giving a ratio of the freestream static' pressure to 
vent chamber static pressure of 1.1 o2. d he 30P- 
dON configuration and the locations on the side-wall 
where the Darcy law was applied are sliown in Fig. 2. 
These venting regions approximate tlu" regions em- 
ployed in (onfig. 2 of t he ex[)eriment. 

For simplicity, the side wall venting pattern in the 
CFD study followed grid lines, which is the reason 
for the non-smooth boundary outline near the front 
of th(‘ flap in Fig. 2. The porous semi-cin'ular region 
in front of the slat was not modeled in the (TT3. A 
symmetry (slip) wall boundary condition was em- 
ployed on the side- wall at and upstream of approxi- 
mately 0.8c in front of the wing; thus, the side-wall 
boundary layer was relatively "fresh at the start of 
the wing in the computations. 

The top and bottom walls of the LI P I wind tun- 
nel were not modeled in the current study. A pre- 
vious study* ' explored the effects of iiu'luding these 
walls in 2-D computations. The chief effect was to 
raise the upper surface pressures on the flap and 
on the aft end of the main at high angles of attack 
(yielding lower lift, not c|uantified in the reference). 
There was also an effect on the ('omputed wake posi- 
tions. The current structured grid extended approx- 
imately 15c above, below, in front of, and behind 
the wing, and far field Rieinann boundary conditions 
were a[)plied there. Also, symmetry was assumed in 
t he wind tunnel, and only half of t he span was mod- 
eled by the grid: one side modeled tln^ tunnel wall, 


and the other side used symmetry boundary condi- 
tions to model the tunnel center plane. 

The grid used by CFL3D was the same one em- 
ployed in earlier 2-D ('FD studies, except that it 
was diipli('ated with 33 planes in the spanwise di- 
rection, with viscous clustering near the side-wall, 
d'here were 4 zones (with 1-to-l point matching at 
interfaces) in the grid, with 135,425 nodes ])er plane, 
or 4.47 million total nodes. The minimum spacing 
near solid walls was between 0.2 - 3.1 x l()“'’c. The 
viscous unstructured grids had a far field extent of 
I Ur and 840. 8(53 nodes (no brackets), and 1.35 mil- 
lion nodes (including brackets). The minimum s|)a( - 
ing near solid walls was 1 x l()~'’r. 

4 Results 

4.1 Experiment 

The lift curves foi* the five configurations in the 
current tunnel ent ry are shown in Fig. 3. 4 hese val- 
ues of lift were obtained by integrating the pressures 
from the taps lo('ated along the center of the iiiodel. 
As expected, ( 'onfig. 1 (no wall venting) yields lower 
levels of lift at all angh\s of attack. However, some- 
what surprisingly, the angle of attax'k at maximum 
lift is the same as most of the wall- venting runs 
(n — 21°). It was originally thought that ('orner 
flow separation would cause an earlier stall in this 
case, ('onfig. 2 also exhibited unexpected results: 
its maximum lift o('('urs at a /u(//u c angle of attac'k, 
n > 23°, than all the other runs. Previously, the 
angle of attack for maximum lift had almost always 
come out to be o = 21°. 4 his result suggests tliat 
side-wall venting treatment can inffueiu'c the physi- 
cal process or mechanism (s) responsible for limiting 
maximum lift in the wind tunnel. 

Fig. 4 shows the effect of different venting levels on 
the lift curv<' near maximum lift lor Config. 2. Two 
separate entries ar(’ represented on this j)lot. I he 
initial entry, which used A/> = 8 psi, is the same 
curve sliown in Fig. 3. In a later entry, the A/> was 
varied at each angle of attack. Results using Ap — 8 
psi are inconsistent with the initial entry, indicative 
of a hysteresis effect: in the initial entry, the lift 

continues to increa.se through n — 23°, whereas in 
the other the maximum lift occurs at o = 22°. In 
both the A/> = 6 ])si and Ap =12 psi <’ases, the lift 
continues to increase through o = 23°. 

From the last two figures, it is already (‘lear that 
the flow field near maximum lift is very sensitive, and 
the angle of attack where maximum lift occurs can 
vary l)v 2° or more, depending on the side wall treat- 
ment (although the change in f'/. is less than PX ). 
41ie "widening" venting pattern on the flap seems 
to be ne(‘essary to achieve' higher angle of attack for 
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^ Ljna.V' It is not known why including the slat strip 
((.’on fig. 3) loses this advantage. One question that 
can be asked is vvliether the variation in maxiniuni 
lift is a function of the three-dimensionality of the 
How. Unfortunately, in the experiment we only have 
surfac(' [)ressures to address this issue. Surface oil 
flow, disc'ussed in greater detail below, showed no 
evidence of .separated wing-wall juncture flow when 
wall venting was present. Off-surface flow visualiza- 
tion was not i>ossible in this test. 

Fig. 5 shows spanwise pressures on the flap up- 
per surfax'e at an unstowed chord wise location of 
,r/c — ().9*2o, as a function of angle of attack for 
(‘onfigs. 2 and 1. The plot shows that as the angle 
of attack increases beyond a = Ifl®. the flow field 
bec'omes less and less two-dimensional. It also ap- 
pears tliat at the highest angles of attack, Config. 4 
exhibits somewhat more three-dimensionality than 
(’onfig. 2. In particular, the pressure variation is less 
"flat ' near the tunnel (‘enter ])lane. Fhis increased 
three-dimensionality is consistent with the fact that 
( ’on fig. 4 is unal)le to achieve as high a Ci jnar as 
( ’on fig. 2. 

The pressures and the mass flow rate in the vent- 
ing chambers were measured in this test. Ten pres- 
sures were tak(Ui at various loc'ations behind each 
wall’s plate to check for variations with position. 
Although not shown, these were found to vary by 
less than 0.2 psi at any given angle of attack. Also, 
the pressure levels varied by only a small amount, 
less thaii 0.8 psi, as the angle of attack was varied 
(generally the venting pressure increased at higher 
angles of attack). The mass flow rates as a function 
of angle of attack and configuration are shown in 
V\g. (). When no leading edge slat venting is present, 
th(‘ mass flow rate is roughly 0.7 Ibui/secc Iiududing 
slat vc'nting increases the mass flow rate to approx- 
imatffly 1.5 lb„i/sec. 

Wh(ui side- wall venting is employed, surface oil 
flow shows no indication of fl-D corner flow struc- 
tures near maximum lift conditions. However, wi th 
solid side walls (no venting), significant deviation 
from 2-1) flow at the walls near the corners (‘ould be 
s(‘en at all angles of attack near maximum lift. An 
enhaiu’ed photograph is shown in Fig, 7 at o = 20^. 
There is t urning of the flow oxer both the main ele- 
ment and the flap away from the side wall, and there 
are several coni|>lex flow patterns indicating 5-D flow 
features on the side wall as well. At liigher angles of 
attack, similar patterns exist, except the sizes and 
I o(’ at ions of t he features vary somewhat . Oil was not 
applied at lower angles of attack. 


4.2 Computations 

4.2.1 Structured Oriel Results 

A compilation of the (‘omputed results are shown 
in Fig. 8 for Config. 1 (no venting) and Fig. 9 for 
Config. 2 (venting). In ea(*h figure, the experimen- 
tal results are shown as solid symbols. 2-D nf‘sults 
using CFL3D are shown for reference as a solid line, 
and 3-D results from the current study using C’FL3D 
are shown as open symbols connected by lines. All 
lift levels from the 3-D computations are obtained 
by integrating the pressures in the tunnel center 
plane. When massive side-wail effects occur in the 
(-FD (past tlie ("omputed CL ^iar)^ fhe flow gener- 
ally goes unsteady. However, because steady-state 
time marching is employed, only representative lift 
coefficient valm^s are shown to indicate that the an- 
gle of attack for maximum lift has I>een excf^eded in 
the (’ompu tat ions. 

As shown in Fig. 8, when there is no venting mod- 
eled in the C’FD, lift levels are lower than tlu' 2-D 
levels. This is because 3-D flow features in the wall- 
juncture region lower the wing circulation. Results 
agree well with experiment up to roughly o — 15^. 
but above this the (^FD predicts massive side- wall 
effects and the lift drops dramatically. C’learly ( ’FD 
is not showing the same charax'ter as t he experiment 
at higher angles of attack. Surface streamlines at 
a = 19^ are shown in Fig. 10. C^ompare with Fig. 7 
(at a slightly higher angle of attack). Although re- 
sults show similar character on the main element, 
flow over the flap and side wall is different. CVD in- 
dicates only small deviations from streamwise flow 
on the flap, whereas experiment show^s a large sep- 
arated region. However, the surfact streamlines in 
this flow field are deceptive. The off-hody stream- 
lines for this case are shown in Figs. 1 1. Therf' is a 
significant region of 3-D flow oc’curring off tlit sur- 
face above and behind the flap, the general shape of 
which shows a similar footprint to the photograph. 
In cases such as this one, for which a large off- body 
vortex is present , the current grid is probably too 
coarse to adequately resolve the flow feature. Fig. 12 
shows vorticity contours along with a view of the 
spanwise grid in a vertical plane above the trailing 
edge of the flap, (irid underresolution likely con- 
tributes to over-spreading of the vortex and early 
stall in the ('omputatioiis. 

When venting is modeled in the C’FD. results 
agree with experiment at lower angles of attat'k, 
but are over predicted compared to the ex])eriment 
at. higher angles of attack (see Fig. 9). (’ompu ted 
(f.mar changes depending on the magnitude of the 
suction. Using a = U.l, o('curs at a = 19^, 

cr — 0.2 gives Cijnar — 2U, and a = 0.3 gives 


4 

American Institute of Aeronauti<'s and Ast ronauti('s 



Table 2. Siimiiiary of (’FI) venting characteristics 
for C'onfigs. 1 and 2 


cr 

mass flow, Ibni/sec 

predicted (T,ma.r< !l 

0.0 

n/a 

^ 10 

0.1 

O.-W 

19 

0.2 

i.o:? 

21 

0.3 

1 ,:{8 

22 


('ijuar — '2^^- The higher the suction level, the 
closer the '.]-D computation mimics 2-D com[)uta- 
tional results at the highest angles of attack. These 
results are summarized in 'ral)le 2. along with the 
corres|:)onding mass flow rate computed to be exit- 
ing through the ])orous wall as a result of the Darcy 
]:>ressiire-velocity law boundary condition. (V)ni[>ar- 
ing these numbers with the mass flow rate from 
the exi)eriment for (’onfig. 2 (of approximately 0.7 
Ibni/sec), it appears that a rr between 0.1 aiid 0.2 in 
the ('FD would best rej)resent the actual test. All 
computed results with venting for the rest of the pa- 
per use rr = 0.2. 

Now let us examine results for several different, 
angles of attax'k. Fig. 13 shows spanwise pressun^ 
('oefficients on the flap at o = 4°, botli without and 
with side-wall venting. In the case of no venting, 
CFD accurately predi('ts the experimental level near 
the tunnel center plane, but the spanwise variation is 
overpredicted. ( ’FD res\dts on a (’oarser grid are also 
shown in the figure (grid using every other point in 
each coordinate direction); they exhibit little differ- 
eiu'e from results on the finer grid. When vetiting is 
simulated, C’FD results again agree with the exper- 
imental results near the tunnel center plane; in this 
case the spanwise variation is small for both CTT) 
and experiment. 

Fig. 14 shows ('hordwise pressure coeffi(‘ients at 
o — 4^ for Config. 1 (no venting). CFD results are 
giv(m at four span locations, and are ('ompared with 
experimental results at the tunnel center plane (50% 
s]:)an). 2-D computed results are also shown for ref- 
erence, CdT) agrees well with the experiment at 50% 
span. The changes with sj)an station exemplify the 
significant effects of the 3-D flow o[i the surface ])res- 
sures near the back of t he main element and on the 
flap. Fig. 15 shows chord wise pressure coefficients 
at n = 4*^ for venting C’onfig. 2. Agreement at 50% 
span between (’FD and experiment is again excel- 
lent, and results in this ca.se also agree with 2-1) 
('omputations. There is no notT’eable spanwise vari- 
ation in the computed results. 

Figs. 10, 17, and 1^ show spanwise and chord- 


wise ju'essure coefficients at o = 10'^, without and 
with venting. In this (’a.se, (TD results are not 
grid-indepenchuit between the medium and fine grids 
(Fig. 10). Also, the actual Cp levels for CFD and ex- 
periment, at ,r/c = 0.025 on the flap do not agree: 
(’FI) yields somewhat lower levels in general. Sim- 
ilar to results at n = 4^, the (’FD results again 
show spanwi.se variat ion over the flap and baxT of t he 
jiiaiii for the no-venting case (f ig. 17), and very lit- 
tle spanwise variation for the venting ('ase (f ig. 18). 
Other than small differences in the suction j:)eaks 
and over the upper surfax'e of tin* flap and back of 
the main, agreement with ex[)erimental Cp levels is 
fairly good overall. 

Figs. 10, 20, and 21 show spanwise and chord wise 
pressure coeffnaents at a 10^, witliout and with 
venting. Again, the actual Cp levels for C’FD and 
experiment at j'/c = 0.925 on tlie flap do not agree; 
CFD yields lower levels in general. In Fig. 20, (TD 
for the no-venting ('ase exhibits significant spanwise 
\ariation. I his case has already stalled, and surfa('e 
pressures do not agree at all with expxu'iment, whi('h 
does not reacdi ( ^nOil a higher angle of attack. 

CFD results with venting in Fig. 21 show no no- 
ticeable s[)anwise variation. Agreement with exj)eri- 
ment is generally good, although ( TD yields slightly 
stronger suction peaks than experiment and most of 
the flap upper surface and the ba<'k of the main is at 
slightly lower pressure. Results at higher angles of 
attack, with venting, show agreement/disagreement 
with experiment similar to Fig. 21. Plots of all three 
elements are not shown, but details can be seen in 
close- uj) views near the trailing edge of the main in 
Figs. 22 and 23. Here, the underprediction of pres- 
sure' levels in this region is ('learly seen. 

d he reason for ( TD’s increasing deviations from 
experimental trends near maximum lift is still un- 
known. Even with wall venting modeled, the lift 
either tends to be too high (as for o > 19^ using 
CT — 0.2), or massive separation ota'urs and (/imar 
is rea(4ied too early (as occurs for lower levels of suc- 
tion, e.g., a = 0.1 in Fig. 9). Some of the deviations 
in upper surface pressure on the flap and back of the 
main may be due to the fa<'t that the top and bot- 
tom walls were not modeled.^' Nonetheless, the use 
of ( TD appears to require a delicate balancing act 
of ])araiiieters to achieve a reasonable repn'sentation 
of the true 3-D physics. 

In spite of these diffi(’ulties and uncertainties, it 
can be useful to use the CFD results (with vent- 
ing) to explore the influeiu'e of the three-dimensional 
character of the flow field o!i Cijnar- In Fig. 22 at 
('L,nuu‘ — 21°. little spanwise variation is seen, but 
in Fig. 23 at n = 22° (past the predicted Cijna.r)^ 
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s])an\vis(' variation is clearly starting to show up. 
rigs. 21, 25, anti 26 show ofi-l)ociy streamline's for 
flow fields computed at. three angles of attack with 
wall venting. At the two lower angles of attack u|) to 
and including t V.,n;a.r = 21^, there is off- body tliree- 
dimensional corner tiow. but it is relatively small 
and its effend is localized. At n — 22^. however, 
the tliree'-dimensionality increases. Fig. 27 shows 
liow significantly the character of the spanwise pres- 
sure changes at o = 22*^. (It also shows that the 
Cp values on the flap for the 2-D and 5-D computa- 
tions. although indistinguishable at tlie scale sliow!i 
in Fig. 18, really ar<' slightly different at and above 
a 1()^ at .r/c = 0.925.) 

In looking at Fig. 27. one might be tempted to 
decdare that the How at o — 22*^ is ''2-D enough" 
on the port ion of t he wing near the center plane of 
th(' tunnel (because the curv(' is relatively Hat at 
that location), but this view ignores the effect of 
the corner How on the overall circulation of the wing 
system. When a significant region of c-orner How 
circulation exists, there is a lower overall circulation 
around the wing, and the lift level decreases. In 
otiu'r words, f Y,ma.r ii^ t his case is limited l>y iFD 
effects, and not by the 2-D mechanism conjectured 
in 5 ing et al.^^ 

(’omparing Fig. 27 with Fig. 5 from the exper- 
iment, it is seen tliat the two are noi behaving 
similarly. Whereas the CFD spanwise variation re- 
mains relatively Hat through a = 21°, the experi- 
ment shows more evidence of flow even as low as 
o 19°. This difference probably t'xplains why the 
experiment givt's lower levels of lift as it approaches 
( 'L.mar- F iifort u natelv, we were not suct'essful in 
simulating the <‘onditions u('('('ssary to achieve this 
same degree of three-dimensionality in the CFD. 

It should be noted that many of the A-D CFD runs 
were de|x'ndent on the initial conditions. For exam- 
ple. if a high angle of attack run was restarted from a 
solution for wliicdi t he How was massively separated, 
results t<uided to rt mam iiighly .separated. On the 
other hand, ifrt'started from a mostly-att ached-How 
solution, a liigh angle of attack run had a greater 
likelihood of remaining attaclu'd. This issue of non- 
uniqueness in the CFD solution for these flows is 
troubling. However, all runs for this study were 
condui'ted in 'steady state' mode (i.e,, non-time- 
accurate time mari'hing was employed). It is possi- 
ble that when significant regions of corner flow se[>- 
aration are present, then time-accurate approaches 
must be pursued in order to better represent the 
physics. 


4.2.2 Unstruc tured (bad Results 

The unstructured grid methodology was primar- 
ily eiiq)loyed for the purpose of investigat ing the ef- 
fects of the mounting brackets on the .solution near 
CY,N/«.r- It was easic'r to make an unstruct ured 5-D 
grid t o include brackets using VC RID as oj)po.sed 
to ('reating a structured grid. Ilow'ever, even the un- 
structured grid generation ])rocess had limitations. 
It proved to be too difficult to create a .sufficiently- 
refined viscous grid with fine normal spacing near 
the wing, brackc'ts, and in t he wake regions - in t he 
time allotted. In the end. tlie viscous grids used were 
somewhat too coarse to adec|uately repre.sc'nt the lift 
levels of this configuration (lift levels were about 8 
12% lower than results on the fine structured grid), 
therefore, this part of the' ( T^D st udy should be re- 
garded qualitatively only. 

A view of the wing including mount ing bracket s is 
shown in Fig. 28. T he unstructured grid on the side 
wall is also shown. There were a total of 8 Ijrack- 
ets in tlie tunnel, four for the slat and four for the 
flap. The near- wall brackets (0.77 in. wide) were 
at 10.7% span, and the near-c^enter brackets (0.09 
in. wide) were at .40.8% span. This figure shows 
half of the model. Idistructured-grid results were 
obtained using both inviscid (Euler) equations (on 
grids with inviscid-type spacing near the wing) and 
visc’ous (Navier-Stokes) equations on the grid shown 
in the figure. The side wall was coinjiuted using 
inviscid boundary conditions in both cases. Result- 
ing lift coefficients, from integration of the pressures 
along tlie ('enter plane, are given in Fig. 29. The 
effects of the brackets on the inviscid solution are 
minimal, but the brackets lower tin' lift levels in iUv 
viscous ('omputations by approximately 2 3% . Most 
of th(' effect is felt by the main element. Fhis result 
is interesting in light of the fact that the structured 
CFD results (no brac'kets) teinl to over[>redict the 
lift levels near Cl. mar fhe experiment. However, 
firm (conclusions are not possible without additional 
study. 

('omf)uted streamlines for the o = 20° case with 
brackets, sliown in Fig. 40, indicate .some deviation 
from 2-D flow even relatively far away from tlie 
bracket span locations. Pressure coeffi('ients as a 
function of span location are shown in Fig. 4>1 near 
the back of the main element. A noticeable span- 
wis(’ variation is .seen on the upj)er surface of the 
main and flap. Although not shown, when no brack- 
('ts are present, there is very little variation. 

5 Summary and CoiK^hisioiis 

In an effort to determine if 4-D effects were t he 
root cause of previous disagreement between CFI) 
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and experiments for the 30P-:i()N near maximum 
lift, a combination experimental/CFD investigation 
was undertaken. Unfortunately, even the 3-D (TD 
showed discrepancies from the experiment, espe- 
cially at high angles of attack. In particular. (‘FD 
tended to either predict lift levels too higli, or else 
stalled too early. This discrepancy may have been 
due to inadequate grid resolution (e.g., no attempt 
was made to grid- resolve the off- body (‘orner-flow 
vortex), part to inadequate pliysic'a] representation 
(e.g., not modeling the top and bottom walls, not 
representing the side- wall boundary layer inflow cor- 
rectly, not including mounting brackets, not repre- 
senting tunnel disturbances and asymmetries), and 
part to inadetpiacies in modeling (e.g., turbulence 
modeling, Darcy's pressure- velocity law boundary 
condition). 

However, details aside, this study answered (or be- 
gan to answer) some questions. I’he angh' of attack 
of maximum lift in the experiment was clearly sen- 
sitive to side- wall venting shape. 1 his fact strongly 
suggests that sonu' flow feature near the side walls 
may be res|)onsible for limiting rat her than 

a purely 2-D mechanism. Tlje 3-D (‘FD re.^sults in- 
dicated that this feature was probably an off-body 
corner-flow vortex that was present, at high angles 
of attark even when side- wall venting was ap])lied. 
It is consistent that the ^^videning" venting [)attern 
in the exj^eriment achieved the highest be- 

('ause the larger venting area over the flap would 
have had more of an influence on an off-body vor- 
tex. 

A summary of major points from the experimental 
study were: 

• A new venting [)attern with a '"widening r<'gion 
of suction over the flap increased the angle of 
attack at maximum lift from o = to at least 
23 " . 

• The appli('ation of different levels of suction 
(between t) and 12 psi difference from the 
freestream tunnel total pressure) had relatively 
small eff'ects (less than 1%) on lift levels near 
maximum lift. 

• Above o = !()", there was a significant amount 
of spanwise pressure variation on t he flap. ev<'u 
with side- wall venting present. 

• When no side- wail venting was employed, lift 
coefficient levels dropped nearly uniformly by 
0.1 -0.2, but th(' angle of attack at maximum 
lift remained at n =21". 


• The mass flow rates t hrough the side wall vent- 
ing wer(' (juant ified for the purposes of t his and 
future 3-D (‘FD efforts. 

A summary of major points from the (‘FD study 
were: 

• The fine grid was probably sufficiently fine for 
use at the lowc'r angles of att.acT, but its ade- 
quacy at higher angles of attack was dubious, 
particularly because of undc'rresolution of the 
wall vortex. 

• (‘FD with no side- wall venting c'ould predict 
the character of the experiment at lower angles 
of attack. However, above o = 16", the (TD 
[)redict.ed massive 3-D corner flow features and 
conseciucMh ly a loss of lift not seen in the ex- 
periment. The 3-D flow fc'atures were mostly 
off- body; .sur/r/re streamlines over the flap were 
different in character from oil flow in the exj)er- 
iment . 

• (‘FD with side- wall venting c’ould predict the 
character of the experiment at lower angles of 
attack u|) through o = 16", but (TD (using 
enough suction to avoid early stall) tended to re- 
tain more spanwise two-dimensionality t han ex- 
periment at higher angles of attack. Tunnel c^en- 
ter plane pressures still agreed fairly well with 
experiment, but integrated lift near CLmax 
generally overpredic'ted by as mucdi as r>%. 

• (TD yielded different angles of attack for 
(T.rna.r depending on the si<le-wall suc'tion level 
applic'd. The higher the suction, the closer thc" 
results of the 3-D simulation mimicTed 2-D lift 
levels. 

• The maximum lift achievable by C’FD ap])eared 
to be directly related to the degree of three- 
dimensionality in the flow field. At and be- 
low C/,,rna.r. the off-bodv ("orner flow features 
were relatively small and there was very lit- 
tle spanwise variation of surface [pressures. Be- 
yond ( t, 77 ja.r. the off-body corner flow features 
were larger and spanwise variation of surface 
pressures increased. This behavior was difler- 
ent in character from experiment, which exhib- 
ited large spanwise variation of surfax'e pres- 
sures p no I' io Ciffiar as well as after. 

• Darcy's law boundary condition required a rr 
parameter loetween 0.1 and 0.2 to achie\^e a sim- 
ilar mass flow rate as in the experiment. 
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• An Investigation using unstructured grids in- 
di(‘ated tliat the lift was reduced by approxi- 
mately 2-l\% near Ci jfuir when brackets were 
iiK'luded in viscous comj)utations. 

Aside from the [joints listed above, some gen- 
fU'al conclusions from this work were: (a) 2-D CFD 
should not be expec'ted to agree with tlu^ nominally 

2- D wind tunnel exj)eriment at high lift conditions 
be(‘ause the experiment lost its 2-D character at high 
angles of attack; and (b) 3-D (TD using the cur- 
rent grids and methodology compared well with ex- 
periment at low angles of attack, but did not ade- 
quately model the character of the wind tunnel flow 
field near maximum lift. To improve this deficieiu'y, 
based on our experience we recommend that future 

3- D (TD efforts for this {‘onfiguration include (in 
order of importance): 

• finer grid resolution in the region of the wall 
vortex, and overall fiinu’ resolution for unstriK - 
tured grids 

• lop and bottom walls and mounting brackets 

• better (‘haracterizat ion of the incoming side- 
wall boundary layer 
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Figure 2. Location of side- wall venting in the ('om- 
j)ii tat ions (('onfig. 2). 



away in desired pattern 


as constant 

gap 


CONFIG. 4 STRIP ABOVE NAM ELEMENT AND FLAP 
CONFIG. 5 SAME AS CONRC. 4 WTTH 

ADDtTIONAL STRIP IN FRONT OF SLAT 


Figure 1. Location of side-wall venting in the exper- 
iment (dimensions in inches). 



Figure Lift curves in the experiment as a function 
of ('onfignration. 
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Figure 4. F^ffeet of clifTereDt \entiiig levels near max- 
imum lift. Config. 2. 



Figure 5. Spaiiwise surfa(*e [)ressure eoeffieients on 
the fla[) U[)per surfac'e, x /c = 0.925. 
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a, deg 

Figure 6. Measured mass flow rates through eaxdi 
side wall. 



Figure 7. Photograph with superimposed sketch of 
oil flow pattern on wall, flap, and aft end of main for 
Config. 1, a = 20^. 
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a, deg 

Figure 8. Lift roelfieients for Config. 1. 
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Figure 9. Lift coefficients for Config. 2. 
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Figure 10. Coiiipiited surfcKt streamlines looking 
upstream, Config. 1 (no venting), o — 19^. 



Figure 11. Computed ojf?W>or/y streamlines looking 
upstream, ('onfig. 1 (no venting), o = 19°. 
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iMgure 12. \brticity rent ours and grid near side wall 
in vert ical plane above trailing edge of flap, Config. I 
(no venting), o = 19^. 


Figure 14. (fliordwise surface pressur<' <oeffi(dentj 
for (’onfig. 1 (no venting), o := 4®. 
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Config. 1 , CFD 





Config. 1 , CFD (med grid) 





Config. 2, CFD, o=0.2 
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Figure 1(). Spanwise surface pressure (’oefficieiits on 
the flap upper surface at xjc — 0.92r), o = 1(3®. 
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Figure 18. Chord wise surface pressure (oeflicients 
for Config. 2, ct — ().‘2, o = 
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Figure 17. ( hordvvise surface pressure ('oeffi('ieiilf; 
for Config, 1 (uo venting), n == 16^. 
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Figure 19. Spanwise surface pressure (‘oefhcients on 
the flap upper surface at x/c — 0.925, o = 19®. 
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p«rcent span 12.57. 257. 37.57 507. (canlerllne 1 


parcant span 12.57. 257. 37.57. 507. (cantarltna ) 


Figure *24. ('Oiiiputed off-body streamlines looking 
downstream, ('Onfig. *2. rr = 0.2, a = 10^. 


Figure *2(). Computed off'-body streamlines lookiiis 
downstream, (‘onfig. 2, (t — 0.*2, o = 2*2^. 



Figure 25. ('Omputed off-body streamlines looking 
downstream, Config. 2, rr = 0.2. n = 2C. 



IdgLire 27. (.'omputed spamvise surface pressure ec 
efficients on the flap upper surface at j'/c = 0.92.^ 
( 'onfig. 2, rr =: 0,2, 
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Figure 28. V iew of the wing surfac'e including brack- 
ets in the unstructured grid. 


Figure .‘^U. Computed off- body streamlines using un- 
structured metliod with (^rackets, looking upstream, 
0 = 20 ". 
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Figure 29. (’omputed lift coefhcionts using unstruc- 
tured method. 


Figure 31. Close-up of ('hordwise surface pressurr 
coefficients near trailing edge of main using unstruc- 
tured method with brackets, o =20". 
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